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Abstract This paper addresses the workspace analysis of the orthoghde, a 3-DOF 
parallel mechanism designed for machining applications. This machine 
features three fixed parallel linear joints which are mounted orthogonally 
and a mobile platform which moves in the Cartesian x-y-z space with 
fixed orientation. The workspace analysis is conducted on the bases of 
prescribed kinetostatic performances. The interesting features of the or- 
thoglide are a regular Cartesian workspace shape, uniform performances 
in all directions and good compactness. Interval analysis based methods 
for computing the dextrous workspace and the largest cube enclosed in 
this workspace are presented. 
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1. Introduction 

Parallel kinematic machines (PKM) are commonly claimed to offer 
several advantages over their serial counterparts, like high structural 
rigidity, high dynamic capacities and high accuracy (Treib and Zirn, 
1998; Wenger et al., 1999). Thus, PKM are interesting alternative de- 
signs for high-speed machining applications. 

This is why parallel kinematic machine-tools attract the interest of 
more and more researchers and companies. Since the first prototype 
presented in 1994 during the IMTS in Chicago by Gidding&Lewis (the 
Variax), many other prototypes have appeared. 

However, the existing PKM suffer from two major drawbacks, namely, 
a complex Cartesian workspace and highly non linear input/output rela- 



tions. For most PKM, the Jacobian matrix which relates the joint rates 
to the output velocities is not constant and not isotropic. Consequently, 
the performances (e.g. maximum speeds, forces accuracy and rigidity) 
vary considerably for different points in the Cartesian workspace and for 
different directions at one given point. This is a serious drawback for 
machining applications (Kim et al., 1997; Treib and Zirn, 1998; Wcnger 
et al., 2001). To be of interest for machining applications, a PKM should 
preserve good workspace properties, that is, regular shape and accept- 
able kinetostatic performances throughout. In milling applications, the 
machining conditions must remain constant along the whole tool path ( 
Rehsteiner et al., 1999; Rehsteiner et al., 1999). In many research papers, 
this criterion is not taking into account in the algorithmic methods used 
for the optimization of the workspace volume (Luh et al., 1996; Merlet, 
1999; Ottaviano, (1991)). 

The orthoglide is a 3-axis PKM with the advantages a classical serial 
PPP machine tool but not its drawbacks. It is an optimized version of 
the Delta mechanism defined by Clavel, (1990). Most industrial 3-axis 
machine-tool have a serial PPP kinematic architecture with orthogonal 
linear joint axes along the x, y and z directions. Thus, the motion of the 
tool in any of these directions is linearly related to the motion of one 
of the three actuated axes. Also, the performances are constant in the 
most part of the Cartesian workspace, which is a parallelepiped. The 
main drawback is inherent to the serial arrangement of the links, namely, 
poor dynamic performances. 

The orthoglide is a PKM with three fixed linear joints mounted or- 
thogonally. The mobile platform is connected to the linear joints by 
three articulated parallelograms and moves in the Cartesian x-y-z space 
with fixed orientation. Its workspace shape is close to a cube whose 
sides are parallel to the planes xy, yz and xz respectively. The interval 
analysis is conducted on the basis of prescribed bounded velocity and 
force transmission factors. Interval analysis based method is used to 
compute dextrous workspace as well as the largest cube enclosed in this 
workspace (Merlet, 2000). 

2. Description of the Orthoglide 

Most existing PKM can be classified into two main families. The 
PKM of the first family have fixed foot points and variable length struts 
and are generally called "hexapods". The second family of PKM has 
been more recently investigated and have variable foot points and fixed 
length struts. PKMs of the second family are more interesting because 
the actuators are fixed and thus the moving masses are lower than in 
the hexapods and tripods. 



The orthoglide studied in this article is a 3-axis translational parallel 
kinematic machine and is belongs to the second family. Figure 1 shows 
the general kinematic architecture of the orthoglide. The orthoglide has 
three parallel PRPaR identical chains (where P, R and Pa stands for 
Prismatic, Revolute and Parallelogram joint, respectively). The actu- 
ated joints are the three orthogonal linear joints. These joints can be 
actuated by means of linear motors or by conventional rotary motors 
with ball screws. The output body is connected to the linear joints 
through a set of three parallelograms of equal lengths L = BjC^, so 
that it can move only in translation. The first linear joint axis is parallel 
to the X-axis, the second one is parallel to the y-axis and the third one 
is parallel to the z-axis. In figure 1, the base points Ai, A2 and A3 are 
fixed on the i*'^ linear axis such that A1A2 = A1A3 = A2A3, Bi is at 
the intersection of the first revolute axis ij and the second revolute axis 
ji of the i*'^ parallelogram, and Cj is at the intersection of the last two 
revolute joints of the i*^ parallelogram. When each BiCi is aligned with 
the linear joint axis AiBi , the orthoglide is in an isotropic configuration 
and the tool center point P is located at the intersection of the three 
linear joint axes. In this configuration, the base points Ai, A2 and ^3 
are equally distant from P. The symmetric design and the simplicity of 
the kinematic chains (all joints have only one degree of freedom, fig. 2) 
should contribute to lower the manufacturing cost of the orthoglide. The 
orthoglide is free of singularities and self-collisions. The workspace has a 
regular, quasi-cubic shape. The input/output equations are simple and 
the velocity transmission factors are equal to one along the x, y and z 
direction at the isotropic configuration, like in a serial PPP machine ( 
Wenger and Chablat, 2001). 

3. Kinematic Equations and Singularity Analysis 

We recall briefly here the kinematics of the Orthoglide (See Wenger 
and Chablat, 2001 for more details). 

3.1 Kinematic Equations 

Let 6i and Pi denote the joint angles of the parallelogram about the 
axes ij and jj, respectively (fig. 2). Let pi, p2, ps denote the linear joint 
variables, pi = A^Bi and L denote the length of the tree legs, BiCi. 
The position vector p of the tool center point P is defined in a reference 
frame (O, x, y, z) centered at the intersection of the three linear joint 
axes (note that the reference frame has been translated in Fig. 1 for 
more legibility). 



Let p be referred to as the vector of actuated joint rates and p as the 
velocity vector of point P: 

P = [Pl P2 Psf, 



P 



[x y 



p can be written in three different ways by traversing the three chains 
-^i BiCiP: 



p = nipi + (6'iii + /3iji) X (ci - bi) 
p = U2P1 + (^212 + P2h) X (c2 - b2) 
P = napa + {O^i^ + /jsjs) x (03 - bs) 



(la) 
(lb) 
(Ic) 



where bj and Cj are the position vectors of the points Bi and Cj, respec- 
tively, and rij is the direction vector of the linear joints, for i 



1,2,3. 





Figure 1. Orthoglide kinematic architecture 



Figure 2. Leg kinematics 



3.2 Singular configurations 

We want to eliminate the two idle joint rates 9i and (3i from Eqs. (la- 
c), which we do upon dot-multiplying Eqs. (la-c) by Cj — bji 
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Equations (2a-c) can now be cast in vector form, namely 

Ap = Bp 

where A and B are the parallel and serial Jacobian matrices, respec- 
tively: 
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with rji = (cj — bi)-^nj for i = 1, 2, 3. 

Parallel singularities (Chablat and Wenger, 1998) occur when the de- 
terminant of the matrix A vanishes, i.e. when det{A) = 0. In such 
configurations, it is possible to move locally the mobile platform whereas 
the actuated joints are locked. These singularities are particularly un- 
desirable because the structure cannot resist any force. Eq. (3a) shows 
that the parallel singularities occur when: 

(ci - bi) = a(c2 - b2) -I- A(c3 - bs) 

that is when the points Bi, Ci, B2, C2, -B3 and C3 are coplanar. A 
particular case occurs when the links -BjCj are parallel: 

(ci-bi)//(c2-b2) and (c2-b2)//(c3-b3) and (c3-b3)//(ci-bi) 

Serial singularities arise when the serial Jacobian matrix B is no longer 

invertibic i.e. when det(B) =0. At a serial singularity a direction exists 
along which any cartesian velocity cannot be produced. Eq. (3b) shows 
that det(B) = when for one leg i, (bj — aj) _L (cj — bj). 
When A and B are not singular, we obtain the relations, 

p = Jp with J = A~^B (4) 

3.3 Velocity transmission factors 

For joint rates belonging to a unit ball, namely, \ \p\\ < 1, the Cartesian 
velocities belong to an ellipsoid such that: 

p^{33^)p < 1 

The eigenvectors of matrix 33^ define the direction of its principal axes 
of this ellipsoid and the square roots ipi, ^2 and of the eigenvalues 
(7i, (72 and (73 of JJ-^, i.e. the lengths of the aforementioned principal 

axes arc the velocity transmission factors in the directions. To limit the 
variations of this factor in the Cartesian workspace, we set 

'4'min ^ V"* ^ '4'max (5) 

throughout the workspace. To simplify the problem, we set V'min = 
l/tpmax where the value of ipmax depends on the performance require- 
ments. 



4. Determination of the dextrous workspace 

The dextrous workspace W is here defined as the loci of the points for 
which all the eigenvalues of the matrix JJ^, i.e. the velocity transmis- 
sion factors, lie within a predefined range [Cfnim ^max 

]. These eigenvalues 



are determined by solving the third degree characteristic polynomial of 
the matrix which is defined only for the points within the intersection X 
of the three cylinders defined by 

+ < L + <L y'^ + z'^ <L (6) 

To solve numerically the above equations, the length of the legs is nor- 
malized, i.e. we set L = 1. 

Our purpose in this section is to determine an approximation of W as 
a set of 3D Cartesian boxes for any point of which we are sure that the 
constraints on the eigenvalues are satisfied. The width of all the boxes 
in the list will be greater than a given threshold and the value of this 
threshold will define the quality of the approximation. 

4.1 Box verification 

A basic tool of the algorithm is a module M.[B) that takes as input 
a box B belonging to T and whose output is: 

■ either that for any point in the box the eigenvalues lie in the range 

i'^minj '^max\ 

■ or that for any point in the box one of the eigenvalues is either 
lower than (Jmin or larger than Gmax 

■ or that the two previous conditions does not hold for all the points 
of the box i.e. that for some points the eigenvalues lie in the range 
['^min-,(^max\ whilc this is not truc for some other points 

The first step of this module consists in considering an arbitrary point 
of the box (e.g. its center) and to compute the eigenvalues at this point: 
either all of them lie in the range [(Jmin, <^max\ or at least one of them lie 
outside this range. 

In the first case if we are able to check that there is no point in B such 
that the eigenvalue at this point may be equal to a„iin or amax, then 
we may guarantee that for any point of B the eigenvalues will be in the 
range [amin, f^max]- Indeed assume that at a given point of B the lowest 
eigenvalue is lower than amin- this implies that somewhere along the 
line joining this point to the center of the box the lowest eigenvalue will 
be exactly (Jmin- To perform this check we substitute in the polynomial 
A successively by (Jmin and Umax to get a polynomial in x,y,z only. We 
have now to verify if there is at least one value for these three variables 
that cancel the polynomial, being understood that these values have to 
define a point belonging to B: this is done by using an interval analysis 
algorithm from the ALIAS library (Merlet, 2000). 



Assume now that at the center of the box the largest eigenvalue is 
greater than amax- If we are to determine that there is no point of B 
such that one of the eigenvalue is equal to amax , then we may guarantee 
that for any point of B the largest eigenvalue will always be greater than 
(^max- This check is performed by using the same method than in the 
previous case. Hence the Ai module will return: 

■ 1: if for all points of B the eigenvalues lie in [cTmim f^max] (hence 
B is in the dextrous workspace) 

■ -1: if for all points of B either the largest eigenvalue is always 
greater than amax or the lowest eigenvalue is lower than amin 
(hence B is outside the dextrous workspace). 

■ 0: in the other cases i.e. parts of B may be either outside or inside 
the dextrous workspace 

4.2 Algorithm for the determination of the 
dextrous workspace 

The principle of the algorithm is pretty simple: we will maintain a 
list £ of boxes, indexed by i, which is is initialized with the box [-1,1], 
[-1,1], [-1,1]. A minimal width e for the ranges in a box is defined and 
the operator W{Bi) will return the largest width of the ranges in Bi. 
An error index £ will be computed as the total volume of the boxes that 
are not in the approximation but may contain points that are inside the 
dextrous workspace. We then apply Ai{Bi): 

■ if M{Bi)=l: wc store Bi as part of the dextrous workspace and 

consider the next box in C 

■ if M.{Bi)=-l: we consider the next box in C 
m ifM{Bi)=0: 

— if VV(-Bj) > e: we create 2 new boxes from Bi by bisecting the 
range of Bi with the largest width. The two new boxes are 
put at the end of C 

— otherwise we add the volume of Bi to £ 

The algorithm stops when all the boxes in C have been processed. Note 
that this basic algorithm has to be modified in order to consider only 
boxes that belongs to I but this can be done using the same principle. 
The algorithm returns a description of the dextrous workspace as a list 
of boxes and the comparison between the volume of the approximation 
and £ allows to determine the quality of the approximation. Note that 



E is very conservative as part of this volume consists in points that do 
not belong to X or to the dextrous workspace. 

4.3 Implementation and results 

The previous algorithms has been implemented in Maple with a sys- 
tem call to a C++ program that implements the M module. For an 
e of 0.05 we found in about 5 hours that the volume of the dextrous 
workspace is 1.468 with an error bound of [0,0.48] with Omin = 0.25 and 
o-max = 4, i.e. i^rnin = 1/2 and ipmax = 2 (Fig. 3). 

5. Determination of the largest cube enclosed in 
the dextrous workspace 

For usual machine tools, the Cartesian workspace is generally given 
as a function of the size of a right-angled parallelepiped. Due to the 
symmetrical architecture of the orthoglide, the Cartesian workspace has 
a fairly regular shape in which it is possible to include a cube whose 
sides are parallel to the planes xy, yz and xz respectively (Fig. 4). 




Figure 3. Dextrous workspace Figure 4- Cartesian workspace and isotropic 
of the Orthoglide mechanism configuration of the Orthoghde mechanism 

We will now describe a method for determining a cube that is enclosed 
in this workspace whose edge length is 2R such that there is no other 
cube enclosed in the workspace with an edge length of 2(i? + a), where 
a is an accuracy threshold fixed in advance. 

The first step is to determine the largest cube enclosed in the workspace 
with a center located at (0,0,0). This is done by using the A4 module 
on the box [—A — ka, A + ka] , [—A — ka, A + ka] , [—A — ka, A + ka] 
where k is an integer initialized to 1 and A a number initialized to 0. 



Each time the Ai module returns 1 (which means that the cube with 
edge length 2{A + ka) is enclosed in the dextrous workspace) we set A 
to A + ka and we double the value of k. When this module returns -1 if 
k is greater than 1 we reset fc to 1 and restart the process, otherwise the 
process stops and we have determined that the cube with edge length 2A 
is enclosed in the workspace while the cube with edge length 2[A + a) 
is not: thus A is an initial value for R. 

We use then the same algorithm than for the determination of the 
workspace with the following modifications: 

■ for each box Bi in £ we test if at the 8 corners of the cube cen- 
tered at the center of the box and with edge length 2[R + a) the 
eigenvalues satisfies the constraints: 

— if this not the case let u be the maximal half-width of the 
ranges of and assume that ui = R — u > 0. If for at 
least for one corner of the cube C with center of B^ and edge 

length ui the eigenvalues do not satisfy the constraint, then 
the center of the largest cube cannot be located in Bf. indeed 
any cube with a center in Bi and edge length 2R will include 
all the corners of C and thus contain at least a point for which 
the eigenvalues do not satisfy the constraints. 

— if at the corners the eigenvalues are all valid then we search 
for the largest cube centered at the center of the box using 
the same method than for the determination of the largest 
cube centered at (0,0,0). If we find a cube with a larger edge 
length than 2R, then the value of R is updated. 

■ the boxes are bisected only if the largest width of their ranges 
is greater than 2a: indeed even if the largest cube has a center 
located in such box, then the maximal edge length will be at most 
2{R + a). 

Using this algorithm with a value of 0.001 for a we found out that the 
largest cube has its center located at (0.085938, 0.085938, 0.085938) and 
that its edge length was L^Y^j^gpace = 0.643950 while we may guarantee 
that there is no cube with edge length larger than 0.643952 enclosed in 
the workspace 

6. Conclusions 

The dextrous Cartesian workspace and the largest cube enclosed in 
this workspace is computed using interval analysis based method for 
the Orthoglide. Unlike most existing PKMs, the dextrous workspace is 
fairly regular and the performances are homogeneous in it. Thus, the 



entire workspace is really available for tool paths. The bounds of ve- 
locity and force transmission factors used in this paper are given as an 
example because it depends on the performance requirements for ma- 
chining applications. A l:3-scale prototype of this mechanism is under 
construction in our laboratory with these bound contraints. 
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